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Magnetic-field asymmetry of electron wave packet transmission in bent channels 

capacitively coupled to a metal gate 
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We study the electron wave packet moving through a bent channeL We demonstrate that the 
packet transmission probabihty becomes an uneven function of the magnetic field when the electron 
packet is capacitively coupled to a metal plate. The coupling occurs through a non-linear potential 
which translates a different kinetics of the transport for opposite magnetic field orientations into a 
different potential felt by the scattered electron. 
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The current (/) that flows through mesoscopic conduc- 
tors in contact with electron reservoirs is expressed as the 
product of the voltage {V) applied to the leads and the 
conductance (G) which in the linear transport regime de- 
pends on the external magnetic field (B) but not on the 
voltage, I{B,V) = G{B)V. According to the Landauer 
theory the linear conductance (G) of a two-terminal con- 
ductor G{B) = j^Ti2{B), is proportional to the transfer 
probability of a Fermi level electron from terminal 1 to 
terminal 2 (Ti2). For a conductor with disorder (or in 
general with a potential landscape) that is independent 
of B, or is an even function of B, the transmission prob- 
ability is symmetric with respect to the magnetic field 
orientation. The relation 



Ti2{B)^Ti2{-B) 



(1) 



known as Onsager symmetry when expressed in terms 
of conductance, can be derived in the scattering ma- 
trix theory [2| . Time- reversal symmetry alone gives only 
Ti2{B) = T2i{-B), which is due to the fact that the 
electron trajectory from terminal 1 to 2 deflected by 
the Lorentz force for magnetic fleld +B coincides with 
the trajectory from 2 to 1 and opposite field orientation 
—B. The symmetry (1) results from the backscatter- 
ing probability being an even function of the magnetic 
field Rii{B) = Rii{—B), since the backscattered trajec- 
tories are independent of the magnetic field orientation 
[see Fig. 1(a)]. Recently it was pointed out [3| that a 
current asymmetry I{B, V) ^ I{—B, V) can occur in the 
non-linear transport regime, i.e., when the conductor is 
driven out of equilibrium. Magnetic-field asymmetry of 
the non-linear conductance was indeed found in carbon 
nanotubes [J] , quantum billiards [5| rings [6( and dots [7| . 
In this Letter we report on simulations of the trans- 
port of an electron packet through a bent semiconductor 
channel that is capacitively coupled to a metal plate. In 
gated semiconductor nanodevices [8| the confined charge 
carriers are capacitively coupled to the electrodes. The 
coupling is described by the Schrodinger-Poisson scheme 
[9|. The electron wave packet induces a redistribution 
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FIG. 1: (color online) (a) Backscattered trajectories for oppo- 
site magnetic field orientations. For —B {+B) the trajectories 
are deflected to the right (left) by the Lorentz force, (b) Snap- 
shots of the wave packet density along the axis of the straight 
wire that is covered (blue lines) and uncovered (red lines) by 
a metal plate, c) Dots show the position of the centers of 
the basis function [Eq. (2)] deflning the channel. Shaded 
area shows the channel range determined as the region where 
the sum of the absolute values of the basis functions is larger 
than 6.5% of its maximum value. The horizontal lines show 
the range of the bent part y £ [330nm, 726nm] . 



of the charge on the surface of the metal. The redis- 
tributed charge on the metal surface is a source of poten- 
tial that focuses the electron wave function and stabilizes 
the shape of the packet when it is set in motion [lO| . The 
electron density interaction with the metal is non-linear 
and results in a cancelation of dispersion effects intrinsi- 
cally present in the linear Schrodinger equation |ll| . For 
that reason the electron wave packet interacting with the 
metal has a soliton character. We demonstrate below 
that the transmission probability of the electron soliton 
through a bent channel is an uneven function of the mag- 
netic field. 

Our simulation is based on a numerical solution of 
the two-dimensional time-dependent Schrodinger equa- 
tion ih^ = H^ . We use an approach 12] originally 
introduced to describe the attenuation of the Aharonov- 
Bohm oscillations by the preferential injection of the in- 



coining wave packet to one of the arms of the quantum 
ring by the Lorentz force. The kinetic energy Hamil- 
tonian is given by T = {—ihV + eA) /2m*. The Lan- 
dau gauge A — {—By, 0, 0), GaAs electron effective mass 
TO* = 0.067too and dielectric constant e = 12.4 are used. 
The time evolution is calculated in the basis of Gaussian 
functions *(r, f) = Z]„ c„(t)/„(r), where 

/„(r) = exp{-a(r - R^f + iP{x - X„)(y + y„)}, (2) 

with r = (a;,y), (3 = eB/2h, a = l/(22nm)2 and 
R„ = {Xn,Yn) is the center of the n-th Gaussian. The 
imaginary term in (2) is related to the magnetic trans- 
lation and ensures gauge invariance. In the applied ap- 
proach the centers R„ define the structure itself. The 
choice modeling the bent channel is shown in Fig. 1(c). 
We assume that at a distance d = 20 nm above the 
channel the structure is covered by an infinite metal 
plane. For infinite plane the solution of the Poisson equa- 
tion [10] for the response of the metal can be obtained 
with the image-charge technique [11]. Interaction energy 
of the electron density with metal introduces is therefore 
given by 



0.75T -0.75T 

i i 



0.75T -0.75T 



Vir,t) 



le*(r',t)l2dr' 



from the Crank- 
At) - *(t)] /At = 
The Hamiltonian 



47reeo [{x - x')^ + {y - y'Y + {2dY]^'^ ' 

(3) 
All the time dependence of the wave function (2) 
is introduced in the expansion coefficients c„(t). 
Equations for c„(i) are obtained 
Nicolson differential scheme ih [^(t 4 
[H{t + At)*(i -I- At) + H{t)^{t)] /2. 
depends on time through the induced potential H{t) = 
T + V{t). For the basis (2) the Crank-Nicolson scheme is 
given by a system of equations Sc(t + At) = Sc(t) — 
^[H(t-h At)c(t-h At)-HH(t)c(t)], where the overlap 
and Hamiltonian matrix elements are given by Smn = 
ifnilfn) and H,„„(t) = {fm\H{t)\fn), respectively. As 
the initial condition we take a self-focused wave packet 



11) localized 500 nm before the bent segment moving 
121 ] along the incoming wire with average kinetic energy 
n'^P/2m* = 1.4 meV (fc = 0.05/nm). 

A comparison of the motion of the soliton and normal 
(i.e. non-interacting with metal) wave packet moving in 
a straight channel is shown in Fig. 1(b). Scattering on 
the bent segment is illustrated in Fig. 2 which shows the 
density of the moving packet at specific times (increasing 
from top to bottom figures) for the situation without 
(two left columns) and with a metal gate present (two 
right columns) for B = ±0.75T. For positive value of 
the field the wave packet is deflected to the left by the 
Lorentz force as it moves up [see also Fig. 1(b)] and 
thus it is injected into the bent of the wire [Figs. 2(a,e)]. 
For negative B value the injection occurs only after the 
velocity of the packet is inverted [Fig. 2(e,m)] when the 
packet is reflected at the entrance. This introduces a 
delay in the transport of the packet for negative B value. 
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FIG. 2: (color online) Snapshots of the motion of the electron 
packet injected into the bent wire from the lower lead for 
opposite magnetic field orientations. Later moments in time 
correspond to lower plots. Solution of the linear Schrodinger 
equation (a-h) is presented for t — 2.1, 4.8, 6.7 and 9.5 ps and 
the solution including the non-linear term of the wave packet- 
metal interaction (i-p) is given for t = 3.4, 5.8, 7.7 and 9.7 ps. 
Green lines delimit the channel as defined in Fig. [ijc). 



More details of the motion are seen on the plots of the 
density taken along the straight part of the wires {x = 0). 
For B = —0.75 T we notice [Fig. 3(a)] a compression of 
the wave packet which occurs when the packet is reflected 
at the entrance and at the exit of the bent part. No com- 
pression is seen for B = 0.75 T, for which the Lorentz 
force guides the packet into the bent and out of it with- 
out reflection at the junctions. For both B = ±0.75 T 
the wave packet transfer probability is nearly 95% and 
no reflected part of the packet can be noticed at this 
scale of the plot. The reflected part is visible for zero 
magnetic field [Fig. 3(b)]. In all cases the transmit- 
ted packet stays more or less compact, but with a shape 
that changes in time. The breathing of the transmitted 
wave packet is due to deviations of its form from the sta- 
ble one that occurs during the transfer through the bent 
segment. For B — +0.75 T the breathing is nearly ab- 
sent, and it is strongest for B = —0.75 T for which wave 
packet compression at the junctions occurs. The results 
for the normal packet are given in Fig. 4(a) where the 
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FIG. 3: (color online) The snapshots of the packet density 
along the axis of the straight channels (x — 0) for the elec- 
tron interacting with metal for B = —0.75 T, and 0.75 T in 
(a), (b) and (c), respectively. In (b) the incoming, reflected 
and transmitted wave packets are marked by i, r and t, re- 
spectively. 
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FIG. 4; (color online) Same as Fig. 3 for the normal wave 
packet (a) and a wave packet interacting with the metal plate 
(b). 



densities are calculated at the same moments in time for 
both field orientations. Both densities in the straight 
wire are exactly the same in the region before the bent 
part for the incoming wave packet as well as for the in- 
terference which occurs when slower (still incoming) and 
faster (already reflected) parts of the wave packet meet 
below the bent. Note that such an interference was not 



observed for the soliton wave packet. Differences in the 
shape of the packet for B = ±0.75 T are clearly visible 
at the entrance to the bent (compression for negative B) 
and in the transmitted wave packet. The reflected part 
of the soliton scattering for B = ±0.4 T is observed in 
Fig. 4(b). Only the plots for the incoming wave packet 
coincide for both B orientations. The snapshots of the 
reflected parts of the wave packet are no longer the same 
when the magnetic field is inverted. We observe here a 
destruction of the soliton which is split into two parts. 
Note, that the reflected part of the scattered soliton [see 
Fig. 4(b)] spreads like a normal wave packet. The image 
charge of this part is not large enough to focus it. 

The time-dependence of the fraction of the wave packet 
inside the bent as well as in the straight wires before and 
after the bent [see Fig. 1(c)] are shown in Fig. 5(a) for 
the normal packet and in Fig. 5(b) for the soliton. In the 
absence of the metal the fraction of the wave packet be- 
fore the bent is exactly the same for both ±0.4 T. Since 
the bent segment does not permanently capture any part 
of the normal packet the fraction of the packet which is 
eventually transmitted to the outgoing wire is the same 
for both B orientations. For the wave packet interacting 
with the metal [see Fig. 5(c)] the reflected part of the 
wave packet before the bent segment is different for ±0.4 
T and a different transmission probability is obtained. 
Soliton destruction is inelastic and associated with de- 
creased absolute value of both electron-metal interaction 
energy and the kinetic energy of the packet. 

The fraction of the wave packet which is transmitted 
into the outgoing wire at the end of the simulation is 
a lower bound to the transmission probability. A small 
fraction of the wave packet still bounces back and forth 
within the bent leaving it in small portions to the incom- 
ing and to the outgoing wires. For that reason the sum 
of the transmitted wave packet and the part which is still 
present in the bent at the end of the simulation is an up- 
per bound to the transmission probability. The bounds 
on T{B) for the wave packet interacting with the metal 
plate are plotted in Fig. 5(c). We conclude that the 
transmission probability is: i) an uneven function of the 
magnetic field, and ii) for i? > 0.1 T the asymmetry in 
the transmission probability for opposite B orientations 
exceeds the evaluation accuracy. 

In the absence of the metal plate the fraction of the 
normal packet in the incoming lead R{B, t) is an even 
function of the magnetic field at any moment in time. 
In general also for the normal packet the transferred 
fraction is different for opposite field orientations, i.e. 
T{B, t) ^ T{—B, t), and a symmetric transfer probability 
is obtained only in the large t limit. The kinetics of the 
transfer of both normal and solitary packets through an 
asymmetric channel is different for opposite field orienta- 
tions. The different kinetics for the wave packet interact- 
ing with the metal is translated into a different potential 
felt by the scattered electron. For that reason the ar- 
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FIG. 5: (color online) Parts of the wave packet before, after and inside the bent part of the wire for the electron density 
uncoupled (a) and coupled (b) to the metal plate. Solid lines show the results for B > and dashed lines for B < 0. In (a) 
the part 'before' the packet for B < is given by the symbols, (c) Lower (solid) and upper (dashed) curves are bounds for the 
wave packet transfer probability for B > (blue lines) and B < (red lines). 



gument of the equivalent backscattered trajectories for 
opposite B orientations [cf. Fig. 1(a)] does not hold 
in the present case. The effective scattering potential is 
different when the magnetic field is inverted and conse- 
quently R{B, t) ^ R{—B, t), also in the large t limit. The 
reason of the magnetic-asymmetry is the asymmetry of 
the inelastic soliton destruction for ±B. 

To further convince ourselves, we studied the case of 
the electron dynamics governed by the Gross-Pitaevski 
equation [13| in which the potential (3) is replaced by 
V{r,t) = -7|vl/(r,t)|2 with 7 > (obtained for l/ri2 
replaced by Dirac 6{ri2) in the integrand of (3)). For 
this potential the Hamiltonian matrix elements can be 
evaluated analytically. We found that also in this case 
the transmission probability is asymmetric in B. 

The presented theory was applied to an empty ballistic 
channel that is capacitively coupled to a metal plate with 
a carefully prepared initial wave packet. The channel is 
not connected to electron reservoirs and thus the model- 
ing does not correspond to an electron gas in equilibrium. 
Although the theory does not allow us to calculate the 
current, it is nevertheless related to the previous work 
on the magnetic-field asymmetry of non-linear conduc- 
tance. In the original idea [3|] it was the dependence of 
the potential landscape on B that resulted in the asym- 
metry. In our model the role of the landscape is played 
by the response potential of the metal which depends on 
B orientation through the time dependence of the image 
charge. In previous work [14| it was pointed out that the 
uneven modification of the potential landscape and the 
resulting magneto-asymmetry of non-linear transport is 
due to the Coulomb interaction. In our case the mov- 
ing electron wave packet interacts with the electrons in 
metal. It was recently explained that the Onsager sym- 
metry is broken for conductors interacting with the en- 
vironment [1511. A second conductor [l^] was used as the 



environment which was driven out of equilibrium by the 
applied bias. In the present work the metal is perturbed 
by the moving wave packet. 

In conclusion, we demonstrated that the transfer prob- 
ability of an electron wave packet through a semi-circular 
channel in interaction with a metal plate, that leads to a 
non-linear potential, is not an even function of the mag- 
netic field. 
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